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Abstract 

We study the branching ratios of rare K and B decays in models with minimal flavour viola- 
tion, using the presently available information from the universal unitarity triangle analysis 
and from the measurements of Br(B — > X s j), Br(B — > X s l + l~) and Br(K + — > tt + vD). 
We find the following upper bounds: Br(K + — > ir + vv) < 11.9 x 1CP 11 , Br{K\^ — > 
7r°i/p) < 4.6 x 10 _u , Br(K L -> ^ + M~)sd < 1.4 X 10~ 9 , Br(B -> X s vv) < 5.2 x 1CT 5 , 
Br(B -> X d vv) < 2.2xlCT 6 , Br(B s -> < 7.4xl(T 9 , Br(B d -> < 2.2xl(T 10 

at 95% probability. We analyze in detail various possible scenarios with positive or neg- 
ative interference of Standard Model and New Physics contributions, and show how an 
improvement of experimental data corresponding to the projected 2010 B factory inte- 
grated luminosities will allow to disentangle and test these different possibilities. Finally 
anticipating that subsequently the leading role in constraining this kind of new physics will 
be taken over by the rare decays K + — > n + vD, — > ifivv and B St d — > that are 

dominated by the Z°-penguin function C, we also present plots for several branching ratios 
as functions of C. We point out an interesting triple correlation between K + — > ir + vD, 
B — » X s j and B — » XJ + l~ present in MFV models. 



1 Introduction 



Recently, great experimental progress has been made in the study of Flavour Changing Neutral 
Current (FCNC) decays, leading not only to an impressive accuracy in the extraction of CKM 
parameters from the Unitarity Triangle (UT) analysis ^H], but also to stringent constraints 
on models with extra sources of flavour and CP violation, although an accidental agreement of 
the UT analysis with the Standard Model (SM) cannot yet be excluded [3J|3]. 

One is then naturally led to consider models with Minimal Flavour Violation (MFV) [3], in 
which flavour and CP violation is governed entirely by the CKM matrix jSJ [7j and the relevant 
operators in effective Hamiltonians for weak decays are the same as in the SM. 

As pointed out in 0, there exists a universal unitarity triangle (UUT) valid in all these 
models, that can be constructed independently of the parameters specific to a given model. 
Moreover, there exist several relations between various branching ratios that allow straightfor- 
ward tests of these models. A review has been given in [5]. 

This formulation of MFV agrees with the one of except for the case of models with 
two Higgs doublets at large tan (3, where also additional operators, strongly suppressed in the 
SM, can contribute significantly and the relations in question are not necessarily satisfied. In 
the present paper, MFV will be defined as in jSHH]. 

As reviewed in |8J, this class of models can be formulated to a very good approximation 
in terms of eleven parameters: four parameters of the CKM matrix and seven values of the 
universal master functions Fi(v) that parametrize the short distance contributions to rare 
decays with v denoting symbolically the parameters of a given MFV model. However, as 
argued in [8J, the new physics contributions to the functions 

S(v), C(v), D'(v), (1.1) 

representing respectively AF = 2 box diagrams, Z°-penguin diagrams and the magnetic photon 
penguin diagrams, are the most relevant ones for phenomenology, with the remaining functions 
producing only minor deviations from the SM in low-energy processes. Several explicit calcu- 
lations within models with MFV confirm this conjecture. We have checked the impact of these 
additional functions on our analysis, and we will comment on it in Section |21 

Now, the existence of a UUT implies that the four CKM parameters can be determined 
independently of the values of the functions in Moreover, only C(v ) and D'(v) enter the 

branching ratios for radiative and rare decays so that constraining their values by (at least) 
two specific branching ratios allows to obtain straightforwardly the ranges for all branching 
ratios within the class of MFV models. Analyses of that type can be found in [51 IU1ITT]. 1 

The unique decay to determine the function D'(v) is B — > Xsj, whereas a number of decays 
such as K + — > 7r + z^P, — > -k q vv, B s d — > B — > X s ^vv and Kl — > ir°l + l~ can be used 

1 An alternative approach is to extract from rare decays the relevant Wilson coefficients I12II13IH1] . However, 
since in MFV models these coefficients have nontrivial correlations among themselves, we find it more transparent 
to express the physical quantities in terms of the functions in eq. (| 1 . 1 pi . 
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to determine C(v). The decays B — ► A" Si( ^ + ^~ depend on both C(v) and -D'(u) and can be 
used together with B — ► X s 7 and iT + — > 7r + ^P to determine C(v). 

Eventually the decays K + — > -k + vu and — > -k^vv, being the theoretically cleanest ones 
|15l I16| . will be used to determine C(v). However, so far only three events of K + — > ir + vv 
have been observed |171 1181 ITU] and no event of — > 7r°^P, with the same comment applying 
to B s d — > fi + fi~, B — > Xs dvv and if^ — ► ir°l + l~. On the other hand the branching ratio 
for B — > X s 7 has been known for some time and the branching ratio for — > X s l + l~ has 
been recently measured by Belle [2UI and BaBar |21| collaborations. The latter, combined with 
K + — > 7r + ^P, provide presently the best estimate of the range for C(i?) within MFV models. 

The main goals of the present paper are 

• to calculate various branching ratios as functions of C(v) within MFV models, 

• to determine the allowed range for C(y) from presently available data, 

• to find the upper bounds for the branching ratios of K + — > 7r + z^P, — > 7r°z/P, B s ^ — > 

— > X s and i^i — > vr°/ + / _ within MFV models as defined here, 

• to assess the impact of future measurements on MFV models. 

Our paper is organized as follows. Section 2 can be considered as a guide to the literature, 
where the formulae for the branching ratios in question can be found. In this Section we also 
give the list of the input parameters. In Section 3 we present our numerical analysis of various 
branching ratios as functions of C{v) and their expectation values and upper bounds. A brief 
summary of our results is given in Section 4. 

2 Basic Formulae 

In the MFV models considered here there are no new complex phases and flavour changing 
transitions are governed by the CKM matrix. Moreover, the only relevant operators are those 
already present in the SM. Consequently, new physics enters only through the Wilson coeffi- 
cients of the SM operators that can receive additional contributions due to the exchange of 
new virtual particles beyond the SM ones. 

Any weak decay amplitude can be then cast in the simple form 

A(Decay) = £ B^cd^<W^), *K«) = ^SM + ^New (real), (2.1) 

i 

where Fi(v) are the master functions of MFV models 8 

S(v), X(y), Y(v), Z(v), E(v), D'(v), E\v) (2.2) 

with v denoting collectively the parameters of a given MFV model. Examples of models in this 
class are the Two Higgs Doublet Model II and the Minimal Supersymmetric Standard Model 
(MSSM) without new sources of flavour violation and for small or moderate tan (5. Also models 
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with one universal extra dimension |22l I23| and the simplest little Higgs models are of MFV 
type |2ij. 

In order to find the functions Fi(v) in (|2.2|) . one first looks at various functions resulting 
from penguin diagrams: C (Z° penguin), D (7 penguin), E (gluon penguin), D' (7-magnetic 
penguin) and E' (chromomagnetic penguin). Subsequently box diagrams have to be considered. 
Here we have the box function S (AF = 2 transitions), as well as the AF = 1 box functions 
B uv and B u relevant for decays with vv and 11 in the final state, respectively. 

While the AF = 2 box function S and the penguin functions E, D' and E' are gauge 
independent, this is not the case for C, D and the AF = 1 box diagram functions B vv and 
B u . In phenomenological applications it is more convenient to work with gauge independent 
functions 

X(v) = C(v)+B uU {v), Y(v) = C(v) + B lI {v), Z(v) = C(v) + \d(v). (2.3) 

We have the following correspondence between the most interesting FCNC processes and 
the master functions in the MFV models jSJ : 

K° - K° -mixing (e K ) S(v) 

Bis ~ B\ ,-mixing (AM Sjd ) S(v) 

K — > nvv, B — > Xd gVU X(v) 

K L fiji, B dtS II Y(v) 

K L -» TT°l + r Y(v), Z(v), E(v) 

e', AS = 1 X(v), Y(v), Z(v), E(v) 

Nonleptonic AB = 1 X(v), Y(v), Z(v), E(v), E'(v) 

B^X sl D'(v),E'(v) 

B -» X s gluon E'(v) 

B XJ+l- Y(v), Z{v), E(v), D'{v), E'{v) 

This table means that the observables like branching ratios, mass differences AM ds in B® — 
-Bj] s -mixing and the CP violation parameters e and e', can all be to a very good approximation 
expressed in terms of the corresponding master functions and the relevant CKM factors. The 
remaining entries in the formulae for these observables are low-energy quantities such as the 
parameters Bi, that can be calculated within the SM and the QCD factors t]qqd describing the 
renormalization group evolution of operators for scales [i < Mw- These factors being universal 
can be calculated, similarly to Bi, in the SM. The remaining, model-specific QCD corrections 
can be absorbed in the functions F{. 

The formulae for the processes listed above in the SM, given in terms of the master functions 
and CKM factors can be found in many papers. The full list using the same notation is given in 
|27| . An update of these formulae with additional references is given in two papers on universal 
extra dimensions |221 123j . where one has to replace Fi(v, 1/R) by Fi(y) to obtain the formulae 
in a general MFV model. In what follows we will use the formulae of [221 123) except that: 
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• We will set the functions 

B vp (v), B l \v), E(v) (2.4) 

to their SM values and we will trade the functions D'(v) and E'(v) for the low-energy 
coefficient Cf ff (^) which enters both b — > 57 and b — > sl + l~ . In this manner the only free 
variables are the functions listed in plus the -D(f) function. As remarked below, 

this latter function has only a minor impact on our analysis. We have also explored the 
possible impact of NP contributions to B vv (v) and B ll (v), as will be discussed at the end 
of Section [3 

• In obtaining Br(K + — > 7r + z^p) we have included the recently calculated long distance con- 
tributions 28. that enhance the branching ratio by roughly 6%. This amounts effectively 
to a charm parameter of P c = 0.43 ± 0.07. 

• We will use the formula for (K L -> vr / + r) C p V from [231130"] . 

• We will use the complete NLO formulae for B — > X s j from [31] . 

• We will use the complete NNLO formulae for B — > X s l + l~ from [321 113] . 



Branching Ratios 


Formula 


Reference 


Parameters 


Br(K + -> ir+vv) 


(4.24) 


22. 


Br(K + — > 7T°e + v), m c 


Br(K L -» n uu) 


(4.27) 


22. 


Br(K + -» 7r°e + ^) 


Br(K L -> /iV~)sD 


(4.32) 


22. 


m c 


Br{K\j -► 7r°Z+r) C pv 


(43) 


30. 


see 30 


Br(5 -» 


(4.29) 


22. 


5r(B -» A C ZP) 


£r(5 -> JQz/z?) 


(4.29) 


22. 


Br(B -> X c /P) 


£r(£ s -» 


(4.30) 


22. 




5r(B d -» 


(4.30) 


22. 





Table 1: Guide to the formulae. See text for explanations. The dependence of all branching 
ratios on CKM parameters and the top quark mass is not explicitly reported. 

In Table ^ we indicate where the formulae in question can be found and which additional 
input parameters are involved in them. In Table [2] we give the numerical values of all the 
parameters involved in the analysis. 

Finally, for the reader's convenience, and in order to show the relative importance of NP con- 
tributions to the processes we consider, we report below numerical formulae for the branching 
ratios in terms of in eq. I|2.1|) . These numerical expressions have been obtained for central 
values of the parameters in Table [2J as functions of AC = C(v) — Csm, ACj S = Cj S — Cj$ M , 
AD = D(v)-D SM , AB lI =B lI (v)-B l i M and AB uD = B up (v)-B^. With the aid of eq. (1231) . 
it is possible to quickly check the impact of NP contributions in any given MFV model. As 
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Parameter 


Value 


Gaussian (<r) 


A 


0.2255 


0.0014 




0.0415 


0.0007 


p 


0.191 


0.046 


tj 


0.353 


0.028 


Fb s 


230 MeV 


30 MeV 


FB d 


189 MeV 


27 MeV 


Br(B -> X C W) 


0.1045 


0.0021 


Br(K+ 7r°e + zv) 


0.0487 


0.0006 


pole 


178.0 GeV 


4.3 GeV 


m b 


4.21 GeV 


0.08 GeV 


m c 


1.3 GeV 


0.1 GeV 


a s (M z ) 


0.119 


0.003 



Table 2: Values of the relevant parameters used in the analysis. 

a first insight, we see that the dependence of Br{B — > X s l + l~) on AD is relatively weak, as 
can be read off from the small prefactors in the formulae below. From eq. Q2.5JI one can also 
check whether the NP contribution to box diagrams in any given model is large enough as to 
modify significantly our results obtained for AB 11 = AB VV = in the next Section. Finally, 
these formulae allow to understand the structure of our numerical results. We have 2 : 

Br(B -» X s l + r, 0.04 < g 2 (GeV) < 1) = 1.16 • 10~ 6 (l + 0.38 (AB lT ) 2 + 0.46 ACf AB fl 

+0.41 ACAB a - 3.47 ACf + 0.56 AB lT + 4.31(ACf ) 2 

+0.19 (AC) 2 + 0.38AC - 0.11 AC^AD) , 

Br(B -> X s l + r,l < g 2 (GeV) < 6) = 1.61 • 10" 6 (l + 1.33 (AB ll f + 1.26ACfAB lT 

+1.43 ACAB 11 - 0.31 ADAB 11 + 2.08 AB 11 + 1.42(AC? ff ) 2 

+0.67 (AC) 2 + 1.36AC - 0.29 ACfAD - 0.18AZ?) , 

Br(B -» X s l + r,UA < q 2 (GeV) < 25) = 3.70 • 10~ 7 (l + 1.18 (AB 11 ) 2 + 0.70 ACfAB 11 + 0.60 ACf 

+1.27 ACAB 11 - 0.27 ADAB 11 + 2.18 AB 11 + 0.21(AC£ ff ) 2 
+0.60 (AC) 2 + 1.24AC - 0.16 ACf AD - 0.24AZ)) , 

Br(B d -► n+fT) = 1.08 • 10~ 10 (l + AB fl + Acf , 

Br(B s -► /i+yir) = 3.76 • 10~ 9 (l + AS /[ + AC) 2 , 
2 Notice that we have discarded terms with coefficients smaller than 0.1 in Br(B — > X s l + l~). 
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= 1.50 • 1(T 6 (1 + 0.65 (AC + AB vD )f , 
= 3.67 • 10~ 5 (1 + 0.65 (AC + AB uD )f , 
= 8.30 • 10" 11 (l + 0.20(AC + AB u9 f + 0.89(AC + A-B"*)) 
= 3.10 • 10 -11 (1 + 0.65(AC + AB uD )) 2 , 
= 8.58 • 10~ 10 (1 + 0.82(AC + AB U )) 2 . (2 

3 Numerical Analysis 

Our numerical analysis consists of three steps: 

1. Extracting CKM parameters using the UUT analysis; 

2. Determining the allowed range for AC and AC| ff from presently available data; 

3. Computing the expectation values of rare decays based on these allowed ranges. 

For the first step, we use the very recent results of the JJTfit collaboration on the UUT analysis 
EQ: 

p = 0.191 ±0.046, fj = 0.353 ± 0.028. (3.1) 

Since the UUT analysis is independent of loop functions, the above results are in particular 
independent of the top quark mass. 

In the second step, to minimize the theoretical input, we have traded D'{v) and E'{v) for 
C| ff , which is the relevant low-energy quantity entering Br(B — > X s j) and Br(B — > X s l + l~). 
Concerning Br(B — > X s j), we compare the theoretical value with the experimental results 
of CLEO [331, Belle [2H and BaBar [33] in the corresponding kinematic ranges, adding a 
conservative 10% flat theoretical error to the theoretical prediction. This error contains both 
the uncertainties due to the cutoff in the photon spectrum jHB] and the ones related to higher 
order effects, which are particularly large since we are omitting here model-specific NLO terms 
for the NP contribution. For Br(B — > X s l + l~), we use the experimental data in the q 2 regions 
0.04 < q 2 (GeV) < 1, 1 < g 2 (GeV) < 6 and 14.4 < q 2 (GeV) < 25 to avoid the theoretical 
uncertainty due to the presence of cc resonances. 

The second and third steps are carried out using the approach of ref. jlUJ: taking C(v), 
C| ff (^f,) and D(v) to have a flat a-priori distribution and using the available experimental 
data and theoretical inputs, we determine the a-posteriori probability density function (p.d.f.) 
for C(v), Cf{nb) and all the rare decays listed in Table |3 Concerning D(v), it plays only a 
marginal role in these decays and therefore it is not well determined by the analysis. We varied 
AD in the conservative range ±4Z?sm- Even this rather large variation has little impact on the 
extraction of the allowed range for C{v). 



Br(B -» X d vv 
Br(B -» X a v9 
Br(K + -> vr+i/p 
Br(K L -> vr°i/p 
Br(K L -> u+yu- 
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x10' 




Figure 1: P.d.f.'s for ACf (top-left), AC (top-right) and AC vs. ACf (bottom). Dark (light) 
areas correspond to the 68% (95% ) probability region. 

In Figure[I]we plot the p.d.f. for AC(v) and ACf ff , that represent F^ ew in ()2.1|) and enter 
eq. (|2.5|) , In Figure |21 we plot the p.d.f. for the branching ratios. The corresponding upper 
bounds at 95% probability are reported in Table |HJ where, for comparison, we also report the 
results obtained within the SM, using the same CKM parameters obtained from the UUT 
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Figure 2: P.d.f. 's for the branching ratios of the rare decays Br(K + — > 7r + i/t>) ; Br(K^ — > ifivv), 
Br(K^ — > /i + /i~)sD; Br(B — > XdgVv), and Br{Bd :S ~~ > -DarA; (light) areas correspond 

to the 68% (95%J probability region. 



analysis. Finally, in Figures |3] and 0] we plot the branching ratios of the rare decays vs. C(v), 
to make the impact of future measurements on the determination of C(y) more transparent. 

Let us now comment on our results. As can be seen from Figure we have two possible 
solutions for ACf , one very close to the SM, and the other corresponding to reversing the 
sign of Cf ff (/ife) (recall that 1S negative in the SM and equal to Cf ff (/x&) —0.33). The 

second solution is disfavoured: it is barely accessible at 68% probability, in accordance with 
the results of [Tl|- This result is easy to understand. In the case of the second solution for 
ACf, the branching ratio Br(B -> X s l + l~) b ecomes larger than the experimental value. The 



S 




AC AC AC 

Figure 3: P.d.f. 's for the branching ratios of the rare decays used to constrain AC and ACf ff 
as a function of these parameters: Br(B — » ^ s 7)-B 7 >l.8GeV (top-left), Br(K + — > tt + vi>) (top- 
right), Br(B -> X s l+l~) i4.4<g2( GeV )<25 (bottom-left), Br(B — > X s l + l )i<q2( GeV )<6 (bottom- 
center), Br{B — » ^ s ^ + ^)o.o4<(j 2 (GeV)<i (bottom-right). Dark (light) areas correspond to the 
68% (95%J probability region. Very light areas correspond to the range obtained without using 
the experimental information. 

full results are: 

ACf = (0.02 ± 0.047) U (0.958 ± 0.002) at 68% probability, 

ACf = [-0.039,0.08] U [0.859,1.031] at 95% probability. (3.2) 

Since we have two separate ranges for AC| ff , in the following we will also present separately 
the results corresponding to the "LOW" or "HI" solution for AC| ff (see Figures and • 

As can be seen in Figure ^ we have two solutions for AC, one close to the SM and the 
other corresponding to reversing the sign of C. We recall that C ~ 0.81 in the SM. The ranges 
obtained are 

AC = (-0.16 ± 0.53) U (-2.15 ± 0.08) at 68% probability, 
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Figure 4: P.d.f.'s for the branching ratios of the rare decays Br(K^ — ► -k^vv), Br^K^ — > 
/i + /x~)sD) Br(B — ► Xd, s vv), and Br{B^ s — > /i + /x~) as a function of AC. Dark (light) ar- 
eas correspond to the 68% (95%j probability region. Very light areas correspond to the range 
obtained without using the experimental information. 

AC = [-1.25, 0.44] U [-2.39, -1.45] at 95% probability. (3.3) 
From the plot of AC vs ACf in Fi 

gure^ it is evident that the situation is different for the 
HI and LOW solutions for ACy . Indeed, the two solutions correspond to the following ranges 
for AC: 

LOW : AC = (-0.03 ± 0.41) U (-2.18 ± 0.02) at 68% probability, 
LOW : AC = [-0.75,0.50] U [-2.49, -1.60] at 95% probability, 
HI : AC = (-0.68 ± 0.58) at 68% probability, 

HI : AC = [-1.98,0.04] at 95% probability. (3.4) 

These results are easy to understand. For the LOW solution the solutions with AC being 
positive and negative are consistent with the data on B — > X s l + l~ . On the other hand for the 
HI solution, AC < is favoured as with AC > the difficulty with a too high Br{B — > X s l + l~) 
becomes more acute. 

For the reader's convenience, we report in Table |1] the values of the X, Y and Z functions 
obtained by summing SM and NP contributions and by applying all the available experimental 
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Branching Ratios 


MFV (95%) 


SM (68%) 


SM (95%) 


exp 


Br(K + — > -k+vv\ x If) 11 


< 11.9 


8.3 ± 1.2 


[6 1 1 91 


l 1 ^- ' -8.9 / li^l 


Br(K L — > 7r°iyp) x 10 11 


< 4.59 


3.08 ± 0.56 


[2.03, 4.26] 


< 5.9 • 10 4 [37j 


Br(K^ -> mV~)sd x 10 9 


< 1.36 


0.87 ±0.13 


[0.63,1.15] 




Br(5 -» X s i/P) x 10 5 


< 5.17 


3.66 ±0.21 


[3.25,4.09] 


< 64 38 


Br(B -> x 10 6 


< 2.17 


1.50 ±0.19 


[1.12,1.91] 




Br(B s -> ,u + /-0 x 10 9 


< 7.42 


3.67 ± 1.01 


[1.91,5.91] 


< 2.7 • 10 2 39 


Br(B d -> x !0 10 


< 2.20 


1.04 ±0.34 


[0.47,1.81] 


< 1.5 • 10 3 39 



Table 3: Upper bounds for rare decays in MFV models at 95% probability, the correspond- 
ing values in the SM (using inputs from the UUT analysis) and the available experimental 
information. See the text for details. 



Function 


MFV (68%) 


MFV (95%) 


X 


[-0.71, -0.55] U [0.86,1.90] 


[-0.86, 0.10] U [0.30, 1.95] 


Y 


[-1.23, -1.06] U [0.33,1.37] 


[-1.38, -0.44] U [-0.24, 1.43] 


Z 


[-1.51, -1.40] U [-0.25, 1.31] 


[-1.74, -1.05] U [-0.92, 1.46] 



Table 4: Values at 68% and 95% probability for the functions X, Y and Z. See the text for 
details. 

constraints. 

The impact of Br(K + — > ir + vv) on the bounds on NP contributions can be seen by com- 
paring Figure ^ with Figure [7[ where Br(K + — * 7r + vv) was not used as a constraint. 3 As 
can be seen from Figure the role of Br(K + — > ir + vi>) is to suppress the solution with 
AC ~ —2, which corresponds to destructive interference with the SM in Br(K + — > tt + vv) and 
in the other rare decays. In this respect, a further improvement of the experimental error on 
Br{K + — * n + vv) will be extremely useful in further reducing the importance of this negative- 
interference solution for AC, which is responsible for the peaks around zero for all the rare 
decays in Figure HJ 

We also note that eliminating AC < by means of K + — > 7r + z/t> would basically also 
eliminate the HI solution for ACf ff . We therefore conclude that finding Br(K + — ► ir + vv) 
larger than the SM value would help in eliminating the positive sign of C| ff . To our knowledge 
this triple correlation between K + — > tt + uT> } B — * X s ^ and B — > X s l + l~ has not been discussed 
in the literature so far. It is very peculiar to MFV and is generally not present in models with 
new flavour violating contributions. 

The upper bound on Br(K + — > ■n + ui>) in Table|3]has been obtained using the experimental 
information on this decay. It corresponds to the following 95% probabilty ranges: 

Br(K + -» n + vu) = [0, 0.17] U [0.24, 1.19] 

3 In order to fully exploit the experimental information on Br(K + — > -k + vv), we use directly the likelihood 
function obtained by deriving the experimental CL |41|. 
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Figure 5: P.d.f. 's for the branching ratios of the rare decays Br(K + — > 7r + i/t>) ; Br(K^ — > ifivv), 
Br{K-L — > /i + /x _ )sD; Br{B — ► X^sVv), and Br{B^ s — > /i + /x~) considering only the LOW 
solution for ACy . Dark (light) areas correspond to the 68% (95%,) probability region. 



(LOW: [0,0.12] U [0.39, 1.26], HI : [0, 0.81]) x 10~ 10 . (3.5) 

If we do not use the experimental result on Br(K + — > 7r + i/t>), we obtain instead: 

Br(K + -» 7T+w) = [0, 0.15] U [0.28, 1.12] x 10~ 10 , (3.6) 

corresponding to an upper bound of 11.2 x 10" 11 at 95% probability. 

We have also analyzed the decays Kl — > 7r°e + e~ and Kl — > 7r°/x + /i~ using the formulae 
of |291 130j . In the models with MFV these decays are dominated by the contribution from the 
indirect CP violation that is basically fixed by the measured values of ex and Kg — * ir°l + l~ . 
The dependence on C{v) enters only in the subdominant direct CP-violating component and the 
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B(K%7t*vv)10 10 B(K l ^tc°W)10 11 BR(K L ^ri)10 ! 




Figure 6: P.d.f. 's for the branching ratios of the rare decays Br(K + — > it + vu), Br(K^ — > ifivv), 
Br(K^ — > ^ + /i _ )sD; Br(B — > X^gVv), and Br{Bd :S ~^ considering only the HI solution 

for ACy . Dark (light) areas correspond to the 68% (95% j probability region. 



interference of indirect and direct CP-violating contributions. We find that Br(KL —* 7r°e + e _ ) 
and Br(KL —* 7r°/i + ^~) can be enhanced with respect to the SM value by at most 8% and 
10%, respectively. In view of theoretical uncertainties in these decays that are larger than these 
enhancements, a clear signal of new physics within the MFV scenario is rather unlikely from 
the present perspective. Therefore we do not show the corresponding p.d.f.s. 
Concerning Br^K^ — ► w uu), its 95% probability ranges are given by 

Br(K L -> 7r°i/i7) = [0,4.59] (LOW : [0,4.83], HI : [0,2.84]) x 10" 11 (3.7) 
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x10' 




Figure 7: P.d.f.'s for ACf (top-left), AC (top-right), AC vs. ACf (bottom-left) and 
Br{K + — > tt + ui>) vs AC (bottom-right) obtained without using Br(K + — > tt + uT>) as a con- 
straint. Dark (light) areas correspond to the 68% (95%) probability region. 

(see Figures |2l GO and |f)|) . In Figure |H] we see explicitly the correlation between the charged 
and neutral Kaon decay modes. We observe a very strong correlation, a peculiarity of models 
with MFV I n particular, a large enhancement of Br(Ki, — > 7r°i/p) characteristic of models 
with new complex phases is not possible |43| . An observation of Br^K^ — ► 7r°z/i>) larger than 
6 • 10~ n would be a clear signal of new complex phases or new flavour changing contributions 
that violate the correlations between B and K decays. 
The 95% probability ranges for Br^K^ — > /i + /^~)sD are 

Sr(ir L -► /i + /x") S D = [0,1.36] (LOW : [0,1.44], HI : [0,0.74]) x 10~ 9 . (3.8) 
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Figure 8: P.d.f. for the branching ratios of the rare decays Br{K\, — > n°vu) vs Br(K + — > 
tt + uD). Dark (light) areas correspond to the 68% (9f>%) probability region. Very light areas 
correspond to the range obtained without using the experimental information. 



As in the previous cases, the HI solution corresponds to a much lower upper bound. 
Let us now consider B decays: 

Br{B -» X 8 vu) = [0, 5.17] (LOW : [0, 1.56] U [1.59, 5.4], HI : [0, 3.22]) x 10" 5 , 

Br{B X d vv) = [0,2.17] (LOW : [0,2.26], HI : [0,1.34]) x 10~ 6 , 

Br{B s -» fifi) = [0,7.42] (LOW : [0,7.91], HI : [0,3.94]) x 10~ 9 , 

Br{B d up) = [0, 2.20] (LOW : [0, 2.37], HI : [0, 1.15]) x 10™ 10 . (3.9) 

The reader may wonder whether other observables could help improving the constraints on 
AC and testing MFV models. In particular, the Forward-Backward asymmetry in B — ► X s l + l~ 
is known to be a very sensitive probe of C| ff and of C 01] . Indeed, the HI and LOW solutions for 
ACj S and corresponding possible values of AC give rise to different profiles of the normalized 
Afb, defined as 

S-i d cos °i d2r dsdZ w } s g n ( cos °i) 



^fb(s) 



J-l fltub " dsd cos 6i 



(3.10) 
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This can be seen explicitly in Figure E3 Therefore, a measurement of Afb(s) at a Super B 
factory will be extremely helpful in distinguishing the various scenarios discussed above |45j . 
On the other hand, concerning the CP asymmetry in B —* X s ~f decays |46| . it turns out that in 
MFV models its value is reduced with respect to the SM, once the constraint on the branching 
ratio is taken into account, so that it is not expected to play a significant role in present and 
future analyses [17| . 

In Figure|S]we show the p.d.f.'s for the branching ratios of rare decays for the LOW solution. 
The corresponding result for the HI solution is given in Figure El Clearly the branching ratios 
of various decays are larger in the case of the LOW solution. 




i>- 

0.1 0.2 0.3 OA 0.5 0.6 0.7 O.fi 0.9 1 

a 



Figure 9: P.d.f. 's for the normalized forward-backward asymmetry in B — > X s l + l~ for the LOW 
solution for ACf with AC > —1 (left), for the LOW solution with AC < —1 (center) and for 
the HI solution for ACy (right). Dark (light) areas correspond to the 68% (95%) probability 
region. 



Before concluding this section, let us make a few steps towards the future and consider a 
realistic scenario for the projected integrated luminosities of Belle and BaBar, plus a 10% mea- 
surement of Br(K + — > tt+vv). For concreteness, let us assume the following 2010 experimental 
data: 



Br{B -» X s Z+r) . 04<g2(GeV)<1 = (1.13 ± 0.25) x 10 
Br(B -> X^ + Z")i<g2( GeV)<6 = 



1.49 ±0.21) x 10~ 6 (Belle) 
1.80 ±0.18) x 10" 6 (BaBar) 



Br(B — > X s l + l )l4.4<g 2 (GeV)<25 



(4.18 ±0.48) x 10~ 7 (Belle) 
(5.00 ±0.93) x 10~ 7 (BaBar) 
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Br(B X sl ) 



(3.51 ±0.16) x 10- 4 (Belle) 
(3.67 ±0.16) x 10" 4 (BaBar incl.) 
(3.29 ±0.16) x 10~ 4 (BaBar semincl. 



(3.11) 



corresponding to an integrated luminosity of 1 ab _1 and 700 fb _1 for Belle and BaBar, respec- 
tively. Additionally a reduction to 5% of the theoretical uncertainty in Br(B — > X s j) thanks 
to the ongoing NNLO computation is assumed jlH]. 4 

We can see the dramatic effect of these improvements in Figures 1101121 5-factory data will 
completely eliminate the non-standard solution for ACj S , while they cannot distinguish the 
two solutions for AC (considering only branching ratio measurements), see Figure lT21 However, 
this ambiguity is perfectly resolved by Br(K + — > 7r + i^), leading to the impressive results in 
Figures [TO] and [H] 
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Figure 10: P.d.f.'sforACf (left) and AC (right) in the future scenario specified by eq. 13. lip . 
Dark (light) areas correspond to the 68% (95%) probability region. 



With so powerful experimental data, one can even think of generalizing our analysis by 
allowing for substantial deviations from the SM in box diagrams. If the size of new physics 
contributions to box diagrams is comparable to the SM ones, the results of our "future" anal- 
ysis would not change sizably. On the other hand, a dramatic modification could occur for 
contributions to box diagrams much larger than the SM ones; however, it is very difficult to 
conceive new-physics models in which this possibility can be realized. 

4 Messages 

The main message of our paper is the following one: 

4 The future results for Br(B — * X s j) are referred to the same kinematic ranges as the present results. 
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Figure 11: P.d.f.'s for the branching ratios of the rare decays Br{K + — > 7r + ^z/), Br(K^ — » 
7r°z^) ; BriK-^ — > /i + ^~)sD; Br(B — > and Br(Bd s — > m i/ie future scenario 

specified by eq. h'3.11\) . Dark (light) areas correspond to the 68% (95%,) probability region. 



The existing constraints coming from i<C + — > 7r + ^i>, i? — > X s 7 and 1? — ► X s / + / _ do not 
allow within the MFV scenario of [3] for substantial departures of the branching ratios for all 
rare K and B decays from the SM estimates. This is evident from Table El 

There are other messages signalled by our analysis. These are: 

• The decays B — > X s ^l + l~ will not offer a precise value for the function C even in the 
presence of precise measurements of their branching ratios, unless the theoretical errors 
in these decays and B — > and the experimental error on the branching ratio of the 
latter decay are reduced substantially. This is clearly seen in Figure El 
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Figure 12: P.d.f's for ACf (top-left), AC (top-right) and Br(K + — > n^vv) vs AC (bottom) 
obtained without using Br(K + — > 7r + z^) as a constraint, in the future scenario specified by 
eq. h°J.ll\) . Dark (light) areas correspond to the 68% (95%) probability region. 



• The situation is considerably better in the case of B^^ — » but as seen in Figure |11 
for a given value of Br(Bd iS — ► there are generally two solutions for AC and C, 
that cannot be disentangled on the basis of these decays alone. 

• The great potential of the decays K + — > ir^vv and Kl —* tt°uu in measuring the function 
C is clearly visible in Figures |3] and 01 with the unique value obtained in the case of 
K + — > 7r + z/P in the full allowed range of C. In the case of Kl — ► 7r°^i/ the two solutions 
are only present for Br(Ki — » n°uu) significant smaller that the SM value. Similar 
comment applies to — ► X Sj( 2Z/za 
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• Assuming that future more precise measurements of the K — > ttvu branching ratios will 
be consistent with the MFV upper bounds presented here, the determination of C through 
these decays will imply much sharper predictions for various branching ratios that could 
confirm or rule out the MFV scenario. In this context the correlations between various 
branching ratios discussed in [Sj will play the crucial role. 

• One of such correlations predicts that the measurement of sin 2(3 and of Br(K + — > 7r + z>P) 
implies only two values of Br{Ki — > -n^vv) in the full class of MFV models that corre- 
spond to two signs of the function X 42 . Figure |H] demonstrates that the solution with 
X < 0, corresponding to the values in the left lower corner, is practically ruled out so 
that a unique prediction for Br(Ki — > -k°vv) can in the future be obtained. 

• A strong violation of any of the 95% probability upper bounds on the branching ratios 
considered here by future measurements will imply a failure of MFV as defined in 
unless an explicit MFV scenario can be found in which the contributions of box diagrams 
are significantly larger than assumed here. Dimensional arguments 03] an d explicit 
calculations indicate that such a possibility is rather remote. 

• If the only violation of the upper bounds in Table [31 occurs in B s —> pT and Bj — > 

it will be most likely due to new operators beyond the SM ones. For example, the 
scalar operators which arise in MFV SUSY models at large tan f3 can enhance Br(B s — > 
up to the present experimental upper bound [51 IT51 E9] . 

• Conversely, a violation of the upper bounds for the other channels in Table |3] would 
signal the presence of new sources of flavour and in particular of CP violation. This can 
be confirmed observing a violation of the correlations between K and B decays discussed 
above. 

• In particular, recalling that in most extensions of the SM the decays K — > irvv are 
governed by the single (V — A) g> (V — A) operator, the violation of the upper bounds 
on at least one of the K — > nvv branching ratios, will either signal the presence of new 
complex weak phases at work or new contributions that violate the correlations between 
the B decays and K decays. 

Assuming that the MFV scenario will survive future tests, the next step will be to identify 
the correct model in this class. Clearly, direct searches at high energy colliders can rule out 
or identify specific extensions of the SM. But also FCNC processes can play an important role 
in this context, provided the theoretical and experimental uncertainties in some of them will 
be sufficiently decreased. In this case, by studying simultaneously several branching ratios 
it should be in principle possible to select the correct MFV models by just identifying the 
pattern of enhancements and suppressions relative to the SM that is specific to a given model. 
If this pattern is independent of the values of the parameters defining the model, no detailed 
quantitative analysis of the enhancements and suppressions is required in order to rule it out. 
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As an example the distinction between the MSSM with MFV and the models with one universal 
extra dimension should be straightforward: 

• In the MSSM with MFV the branching ratios for K + -» n+i/p, K L -» ir°vv, B -> X d vv 
and Bd — > are generally suppressed relative to the SM expectations, while those 
governed by like -B — > X s vi?, B s — ► n + fi~ and S — > X^ can be enhanced or suppressed 
depending on the values of parameters involved |5U| . 

• In the model with one universal extra dimension analyzed in |22U23| . branching ratios for 
essentially all rare decays are enhanced, the enhancement being stronger for the decays 
governed by Vt s than for those where Vtd is involved. A prominent exception is the 
suppression of B — ► Xs^j |2*3l IST] . 

Finally, if MFV will be confirmed, and some new particles will be observed, the rare processes 
discussed in this work will constitute a most powerful tool to probe the spectrum of the NP 
model, which might not be entirely accessible via direct studies at the LHC. 
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